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The possibility to induce demixing in a colloidal mixture by adding small polymers, or other 
equivalent depletant agents, is theoretically investigated. By use of Mean Field Theory, suitably 
generalized to deal with short range effective interactions, the phase diagram of a binary mixture of 
colloidal particles (modelled as hard spheres) in a solvent is determined as a function of the polymer 
concentration on the basis of the Asakura-Oosawa model. The topology of the phase diagram 
changes when the relative size of the colloidal particles is reduced: the critical line connecting the 
liquid-vapour critical points of the two pure fluids breaks and the branch starting from the critical 
point of the bigger particles bends to higher volume fractions, where concentration fluctuations 
drive the transition. The effects of a softer colloid-polymer interaction is also investigated: Even the 
presence of a small repulsive tail in the potential gives rise to a significant lowering of the stability 
threshold. In this case, phase transitions may take place by adding just a few percent of depletant 
in volume fraction. These results may be relevant for the interpretation of recent experiments of 
solidification kinetics in colloidal mixtures. 


I. INTRODUCTION 

The realm of liquid state theory [l[, initially developed to study the so called simple fluids like those whose 
constituents are noble gases or small molecules, has been enormously extended when it was realized that such methods 
could be imported to treat many phenomena in a significant fraction of complex liquids, like polymers in solution or 
colloidal suspensions. The strategy is first to identify, among the plethora of the microscopic degrees of freedom, a 
few effective ones and then to construct an effective interaction among these. Models like the hard sphere system or 
the Yukawa fluid found accurate experimental implementations 0 in such context and the analog of the standard 
liquid-vapor and solidification phase transitions are quite common in such complex fluids. However, it should be 
kept in mind that the effective interaction can be very different from those typical of a simple fluid and the phase 
diagram can display more complex topologies compared to the textbook paradigm of three phases, vapor, liquid and 
solid, the last one occupying the high-density side. The liquid phase can disappear altogether or be present only as 
a metastable state Q; re-entrant solidification Q may appear when a fluid region is present at higher density, and 
density modulated phases can also be found Q in place of a liquid or of a solid. 

In the field of simple fluids when one moves from a one component system to a binary mixture a very rich phase 
diagram is found @ and some of them are sketched in Fig. [1] The simplest case is found when the two components 



FIG. 1: Mean field critical lines projected in the plane total packing fraction/volume fraction of component 1, (77 = r;i -|- 772 = 
■|pi-l-f'P 2 , X = for three Lennard-Jones (LJ) mixtures. The diameter ratios 02 ! ai are 0.9 for type I and type II systems 
and 0.8 for the type III model. In the three cases the ratios between the amplitudes of the LJ potentials are: Type I, en = 1.7, 
£12 = 1.3; Type II, en = 1.7, £12 = 1.2; Type III, £11 = 3, £12 = 1.7 in units of £ 22 . The direction of the arrows identify the 
direction of the order parameter: an arrow parallel to the density axis indicates a liquid-vapour transition, an arrow along the 
volume fraction axis, a demixing transition. 
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are rather similar like argon and krypton Q. In a mixture of such atoms one find that the two critical points of the 
pure systems are connected by a line of critical points, well separated in density from the solidification transition; in 
this case, the two components are fully miscible in all the fluid phase. This kind of transition will be referred to as 
“liquid-vapour”. In many, less simple, molecular fluids, like water/phenol mixtures Q, a new phase transition may 
well arise in the high density region, leading to the “demixing” between two liquid phases that differs in concentration. 
This transition ends in a new line of critical points and the phase diagram displays two lines of critical points, one 
in the small density region and one at large density. When the two components are more strongly dissimilar, like 
with Ne-Kr or He-Xe mixtures Q, again one finds two lines of critical points but the topology of the phase diagram 
is different from the previous one because the two critical points of the pure components are no longer connected 
and two disjoint lines depart from them. One line moves toward large density transforming a liquid-vapor transition 
into a kind of demixing one while the other line terminates in an end point. The different character of the order 
parameter at any point along a critical line can be fruitfully identified by an arrow in the density/concentration plane, 
whose direction provides the exact blending of density and concentration fluctuations in the order parameter of the 
transition. Even more complex phase diagrams can be found, and up to six principal types of phase diagrams @ 
have been classified for binary mixtures of small molecules. Not much attention has been given to which extent such 
richness of phase diagrams can be found in the area of colloidal suspensions or of other soft matter systems. This is 
the topic we address in the present paper. 

We consider a mixture of small non-adsorbing polymers and of spherical colloids in a solvent such that the colloids 
can be assimilated to hard spheres. Such polymers give rise to a short-range attraction between the colloids via the 
so-called depletion interaction 0 . When colloids are monodisperse such an attraction gives rise to a phase transition 
between two phases, one rich and one poor in colloids and this transition mimics the liquid-vapor transition of a simple 
fluid, with the concentration of polymers playing the role of inverse temperature. When the range of the effective 
attraction is small enough compared to the colloid diameter, this transition can be metastable with respect to a direct 
transition between a low-density phase and a solid phase at high density of colloids Q. We consider now a similar 
system when colloids are not monodisperse but are a mixture of spheres of two diameters, cti and cr 2 . Since colloids are 
never monodisperse, our model is appropriate when the size distribution of colloids is bimodal with maxima around 
two sizes with a dispersion that is small compared to the difference |cri — (T 2 I. In essence we have a three component 
system, colloids of two comparable sizes and small polymers in a background good solvent. When the polymer degrees 
of freedom are averaged out three depletion interactions arises, the two between colloids of the same size and that 
between colloids of different size. The depletion interaction between two colloids of different diameter is known Q 
in the framework of the Asakura-Oosawa (AO) model in which the inter-polymer interaction is neglected and the 
polymer-colloid interaction is assumed to be hard sphere. However, how such binary AO model affects the phase 
diagram of the system has not been studied and this is our main scope in this paper. In addition, recently, it has 
been pointed out that, by relaxing the hard core hypothesis, a strong enhancement in the strength of the depletion 
interactions can be achieved We have then extended the colloidal mixture/polymer model when polymer-colloid 
interaction is represented by a repulsion over a finite range and it is not infinitely sharp as in the AO model. Based on 
such interaction models we ask if the mixture of two colloids of comparable size has only a liquid-vapor transition like 
in the monodisperse system or if also a demixing transition is present and which is the topology of the phase diagram. 
To this end we have studied the model on the basis of a suitable mean-field approximation and of the mapping onto 
the Adhesive Hard Sphere mixture model 0 determining the phase diagram as function of the size ratio (72/^i, of the 
size of the polymers ^ and of the softness of the polymer-colloid interaction. We find that a demixing phase transition 
can indeed be present and all three topologies of phase diagram displayed in Fig. [T]can be realized depending on the 
model parameters. 

In Section H we first sketch the derivation of the two body effective interaction between two particle immersed 
in an ideal depletant for arbitrary colloid-depletant interaction, thereby generalizing the Asakura-Oosawa model. In 
Section HI we motivate the use of a modified Mean Field free energy to investigate the phase diagram of the colloidal 
suspension. In Section IV the results are shown examining the extension to binary fluids of the Asakura-Oosawa 
model first, and then studying the effects of softer colloid-polymer interactions. Some conclusion and perspective is 
offered in Section V. 


II. THE SOFT POLYMER MODEL: EFFECTIVE INTERACTIONS 

We provide a short derivation of the effective interaction between two colloidal particles, represented by hard-spheres 
of unequal diameters cri and (T 2 respectively, immersed into a non-adsorbing polymer fluid. Following a well established 
procedure 0, polymers are regarded as non interacting point particles, which however feel a repulsive force due to 
the presence of the colloidal particles. Usually such an interaction is represented by the hard core constraint but, in 
view of future generalizations, here we allow for an arbitrary pair potential between a colloid of type i = 1,2 
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and the depletant particles (polymers). The effective colloid-colloid interaction induced by the presence of 

the depletant can be evaluated exactly in this model. The calculation is most easily performed in the semi-grand 
canonical ensemble where the partition function of two colloids, of species i and j, in a reservoir of depletant 
particles is written as: 


y = / dR, dRj Sy(R, - Rj) 


( 1 ) 


where r^ff(R) is the bare potential between the two colloids, (3 = [ksT) ^ and 


N r 


N 


( 2 ) 


is the depletant partition function, with z = ^ the fugacity of the ideal gas, /r its chemical potential and A the 

De Broglie thermal wavelength. The partition function (I2|) is readily evaluated with the result 


\og^ij{Ri -Rj) = z 




( 3 ) 


leading to an effective colloid-colloid interaction 

|3v^J^{R^ - Rj) = /3<;(R* - R,) dr _ i| ( 4 ) 

This result is exact irrespective of the specific form of the bare colloid-colloid (vfj) or colloid-depletant interaction 
v‘[j. Is it customary to write Eq. Q in terms of the “depletant reservoir density” pd which in fact just coincides 
with the fugacity z [ij. By subtracting an irrelevant additive constant to the interaction, the effective inter-colloid 


potential can be written in the final form 


/3</^(R. - R,) = /3C(R. -R^)-Pd dr 


Hr-R.) _ I 


Hr-Rj) _ I 


( 5 ) 


The state dependent additive contribution, omitted in Eq. does not affect the phase behaviour of the mixture, as 
explained in Sec III of Ref. [l^ . We remark that the form of the effective interaction derived here is correct when just 
two colloids are present in the system. In general, many body forces may appear at finite density of colloids [ll|, but, 
when the range of the colloid-depletant interaction is sufficiently small, these terms play a minor role in determining 
the phase diagram of the system and can be safely neglected. The expression for the effective interaction can be 
evaluated analytically in the case of a binary AO mixture, where is a hard sphere interaction with range 

Di = and ^ is a parameter identifying the polymer size, as the polymer gyration radius: 

{ -|-oo for r < Oij 

-pd Vij{r) for CTy <r < (Jij + ^ ( 6 ) 

0 for r > (Tjj -I- ^ 

where aij = ^ 4^0 function vij (r) is given by Q : 


Vijir) 



Q{Dl + D])r + S{Dl + D^^) 



( 7 ) 


In soft matter physics, when the effective potential between colloids is short ranged, a useful parameter measuring the 
strength of the interaction is given by the second virial coefficient or, equivalently, by the so called Baxter’s stickiness 
parameter r 0 defined, in terms of the pair potential v(r) by 


[ dr 

^-Pv{r) _ 

47r o 
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where cr is the hard core diameter. According to the empirical Noro-Erenkel law of corresponding states, the liquid- 
vapour coexistence curve in a fluid interacting via a short range potential is remarkably independent of the details of 
the interaction if the temperature is measured in terms of r. 
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FIG. 2: Demixing ratio as a function of the colloid size ratio at a depletant volume fraction r^d = ~ 0.1. Left panel 

^ = 0.1 (Ti. Right panel ^ = 0.05 ai. Full black curves: AO potential. Dashed blue curves: power law colloid-depletant 
interaction with softness exponent n = 36. Dotted red curves: same with n = 12. 


An estimate of the propensity to demixing of the AO mixture can be given in terms of the stickiness parameters 
of the three effective potentials, as given by Eq. ([8]); If the spatially integrated strength of the unlike interaction is 
smaller then the average of the interactions between particles of the same species, we expect that, at low temperatures, 
a large portion of the phase diagram will display coexistence between phases of considerably different concentration. 
It is then natural to define the demixing ratio in terms of the integrated strengths of the Boltzmann factors of the 
depletion interactions: 


X = 


2 a 


12'12 


. 3.^-1 


+ aKT. 


3.^-1 


( 9 ) 


In the following, the case x < 1, i.e. when demixing is favored, will be named “sub-additive”. A plot of this quantity 
as a function of the size ratio is shown in Fig. [2]for two choices polymer size By lowering the size of the depletant, 
the tendency towards demixing is enhanced. 

We may also consider more general colloid-depletant interactions by relaxing the hard sphere constraint and intro¬ 
ducing a suitable control parameter. As an example, we have investigated a soft core power-law potential parameterized 
by the characteristic length-scales Di for each colloid species i = 1,2 and by the softness exponent n Q: 


(3vfir) = 



( 10 ) 


The colloid-colloid bare interaction is modeled as that of additive hard spheres, leading to the effective potentials 
shown in Fig. [3] for a few representative choices of the parameters. 

The presence of a power-law repulsive decay in the polymer-colloid interaction, even with a remarkably large value 
of the softness exponent n, gives rise to an attractive tail in that may lead to significant changes in the phase 

diagram of the colloidal mixture. To quantify this effect, it is instructive to plot Baxter’s stickiness parameter r as a 
function of the physical depletant volume fraction rjd = for a few choices of the softness exponent. According to 

the analytical solution of the PY equation 0 , a one component fluid is known to undergo a liquid-vapour transition 
when r reaches its critical value Tc ^ 0.0976. In Fig. |4]we can see that the amount of depletant needed to trigger 
the liquid-vapour transition is strongly reduced in going from the original AO model to a colloid-depletant potential 
of the form cni) with n as large as n = 36. A reduction in the polymer size parameter ^ lowers the critical depletant 
volume fraction rid in a roughly linear way. 

On this basis we expect that, also in a colloidal mixture, the occurrence of colloid-depletant interactions of finite 
strength is going to significantly affect the topology of the phase diagram, which is known to be very sensitive to 
the relative strengths of the interactions, as shown in Fig. [21 By increasing the softness of the colloid-depletant 
interaction, the tendency towars demixing is enhanced. 
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FIG. 3: Effective interaction between two hard spheres of diameter ai immersed in a fluid of non interacting depletant particles 
of density p^. The colloid polymer interaction is defined by ^ = 0.1 cti and by a softness exponent n. Full black curve AO 
potential (i.e. n —> oo); dashed blue curve n = 36; dotted red curve n = 12. In the inset the bare colloid-depletant potential 
with n = 36 (dashed blue) is compared with the choice n = 12 (dotted red) 
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FIG. 4: Stickiness parameter r as a function of the depletant volume fraction for three choices of the softness exponent: n = oo 
corresponds to the Asakura-Oosawa interaction, n = 36 and n = 12. The depletant size parameter is ^ = 0.1 cti. The dashed 
line sets the onset of the liquid-vapour transition and the crossing point sets the amount of depletant needed to drive phase 
transition in the colloidal suspension. 


III. FREE ENERGY MODELS FOR A COLLOIDAL MIXTURE 

The theoretical investigation of the phase diagram of binary fluids must rely on some kind of approximation. 
For effective interactions whose range is short with respect to the colloid diameter, standard mean field approaches 
are affected by serious inaccuracies. Even in the simpler case of one component fluids, van der Waals-like theories 
imply the law of corresponding states and lead to quantitative discrepancies with respect to the results of numerical 
simulations Q. The situation is considerably worse when the standard mean field theory is extended to binary fluids, 
because the different topologies of the phase diagram which a mixture may display depend on a delicate balance 
between entropic and energetic contributions to the free energy. Unfortunately, such a balance is not captured by any 
free energy linear in the interactions. In order to circumvent this problem, a simple generalization of the standard 
mean field approach has been rediscovered several times in liquid state theory 0: it simply amounts to assume that, 
outside the core, the direct correlation function coincides with the Mayer function, rather than with the potential 
divided by —fcsT. This simple way to modify the usual Random Phase Approximation [I| (RPA) at least displays the 
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correct low density limit and satisfies the Noro-Frenkel criterion [l^ providing a clear improvement on standard RPA 
with a modest computational effort. By comparing the structure factor predicted by this modified RPA with computer 
simulations or with the results of iutegral equations, a good agreement is found at low and moderate densities, while 
in the high density/low temperature region the modified RPA tends to overestimate the peack height in S{k). 

A different point of view, for dealing with short range potentials, is to resort to the celebrated integral equations of 
liquid state theory, like HNC, PY or MSA [l|. Few of them also allow for an analytic solution in the “adhesive sphere” 
limit, i.e. in the limit of vanishing range of the attractive potential [l^lTsjl. Unfortunately, the analytic solution shows 
that PY does not lead to a full spinodal curve: this integral equation has no solution on one side of the critical point, 
leading to difficulties in extracting the thermodynamics from the knowledge of the correlation function. On the other 
hand, when applied to mixtures, MSA for sticky potentials is known to be affected by a serious inconsistency which 
does not allow to obtain a free energy function via the compressibility route [T^ . Finally, a differential approach to 
the phase equilibria in mixtures (HRT), which proved successful for molecular fluids, has not yet been implemented 
for short range potentials (T^ . 

The phase diagram of our colloidal mixture, interacting via the effective pairwise potential has been 

investigated starting from a mean-field like free energy density of the form: 


^(Pl,P2,Pd) = Ahs{P1tP2) 


keT 
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piPj 


dr 


-/3</UU 



( 11 ) 


where pi and p 2 are the number densities of the two species and Ansipi, P 2 )Js the free energy density of a hard sphere 
mixture, well represented by the Mansoori-Carnahan-Starling expression [l8|. Here we assumed a bare inter-colloid 
potential vfj{r) of the hard core form. This Ansatz, extends in a straightforward way the previously mentioned 
modified RPA for one component fluids [l^. It can be seen as a generalization of the standard Mean Field expression 
for the Helmholtz free energy density: 


AmF {Pl ,P2,Pd) = AhS {Pl, P2 ) 
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piPj 


dr 


„e// 


(r) 



( 12 ) 


where the entropic contribution coincides with that of a pure hard sphere binary fluid, and equals Ahs^ while the 
internal energy is evaluated in Hartree approximation [l|. However, this Mean Field approximation, being linear in 
the attractive part of the interaction can be justified only in the limit of weak and long ranged potential (the 

so called Kac limit), which is not the case of interest (see Fig. |T|). To improve this approximation, maintaining the 
physical content of Mean Field Theory, we can look at the low density limit, where the free energy is known to behave 
as 


A{pi,P2, Pd) —>■ A^'^{pi,p2) 


ksT 
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2 

y^ PiPj 

*.i=i 


dr 



(13) 


The easiest way to enforce this asymptotic behaviour into an approximation of the form (1121) is given by Eq. CD, 
which is in fact our stra tegy . An intriguing direct consequence of such an approximation is that the Noro-Frenkel law 
of corresponding states is here extended to the case of binary mixtures. In fact, within this Mean Field theory, 
the phase diagram is independent of the details of the interactions, provided the temperature is expressed in terms 
of Baxter’s parameters ([5]). Note that when, according to our definition (d), x = 1 the attractive interaction just 
couples the total density pi -|- p 2 in Eq. CD and demixing will not be favoured. 

In the same spirit of Eq. CD, we can evaluate pair correlations in the mixture by use of a modified RPA, expressing 
the direct correlation functions of the mixture as: 

Cij(r) = - 1 (14) 

Starting from the free energy CD, the critical lines are determined by the vanishing of the Hessian determinant of 
A{pi,P2)-- 


det 


d'^A{pi,p2,pd) 

dpzdpj 


(15) 


implying that one eigenvalue of the Hessian matrix vanishes. Let us call and ^2 the two linear combinations of pi 
and p 2 which diagonalize the Hessian matrix at a given point of the spinodal line and let us denote by Ai the vanishing 
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FIG. 5: Projection of the critical lines in the volume fraction/concentration plane for a hard sphere colloidal mixture with an 
AO depletant of diameter ^ = O.lcri. The arrows show the character of the order parameter: an horizontal arrow denotes a 
density fluctuation and a vertical arrow a concentration fluctuation. 


eigenvalue. The critical point is identified by the condition ^ = Oj while the convexity requirement of the 

free energy at the critical point gives rise to two additional conditions: 


' d^A{C,iX2,Pd) 
_ 5Cf5C2 


d'^A{C,iX2,Pd) 

5Cl 

1 gM(Ci,C2,j0d) d‘^A{C,i,C,2,Pd) 
3 SCI SCf 


(16) 

(17) 


These inequalities allow to detect spurious singularities selecting only the stable transition lines. The two “proper” 
variables Ci and C 2 respectively identify the strong and weak fluctuation directions. It is instructive to represent the 
strong fluctuation direction by an arrow in the {pi,p 2 ) plane (T^ . 


IV. PHASE DIAGRAMS OF THE COLLOIDAL MIXTURE 

We numerically investigated the solution of the mean field equations discussed in Section Hill bv varying the model 
parameters. The critical lines and the strong fluctuation direction are shown in the volume fraction vs. concentration 
plane (77 = f (picr? + P 2 ^l)] x = m/v)- 

We first consider a binary fluid of hard particles of diameters cti and <72 with standard AO effective interaction. 
As previously noted, in the pure system (i.e. at concentration a; = 0 or a; = 1) we expect a liquid-vapour transition 
at colloid volume fraction 77 ~ 0.12 and depletant packing fraction rjd larger than rjd ~ 0.26. The projections of the 
critical lines are shown in Fig. [5]for ^ = 0.1 cti and few choices of the colloid size ratio. For a size ratio lower than 0.65 
the transition changes character along the critical line, from liquid-vapour at high concentration of large particles to 
demixing at equimolar concentrations: The phase diagram falls in the class III previously discussed. By increasing 
o’ 2 /o'i, the topology of the critical lines is modified, moving through type II for size ratio 0.7 and finally to type I, 
i.e. pure liquid-vapour, above a 172 0.8 cri. Note that the transition between type II and type I is not sharp and 

the demixing line appears to grow from the high density region, emerging from the solid phase (not shown in these 
diagrams). 

The possibility to obtain an efficient fractionalization of colloids of similar sizes by adding polymers to the suspension 
is however hampered by the large depletant volume fraction required to trigger phase separation {r/d ~ 0.3). We may 
also add that, at those densities, considering polymers as an ideal gas is likely to be rather inaccurate. 

By reducing the polymer gyration ratio to = 0.05 cri, the depletion mechanism is more efficient and the demixing 
transition between the two topologies of the phase diagram occurs at larger colloid size ratio, as can be seen in Fig. [5] 
. Also the amount of depletant required to trigger the instability is lowered down to rjd 0.2, which still corresponds 
to a very dense suspension {rj + rjd 0.4 — 0.5). Reducing the hardness of the colloid-polymer interaction leads to 
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FIG. 6: Same as Fig. [5]for ^ = 0.05 ai. 



FIG. 7: Same as Fig. [5] for a colloid-depletant interaction with softness exponent n = 36 and ^ = 0.1 cri. 


considerable effects on the phase diagram: the sequence of topologies is shifted towards higher size ratios, showing 
that softness enhances the propensity towards demixing, and the required polymer density is reduced, as expected on 
the basis of Fig. |4l As an example, we show in Fig. [7] the sequence of phase diagrams which are obtained for n = 36 
and ^ = 0.1(71. The critical depletant volume fraction along the transition lines varies smoothly around rjd ~ 0.2. Also 
in this case, by reducing the polymer size the depletant density further lowers and the demixing transition occurs 
at a size ratio as high as 0.85, as shown in Fig. [S] Here, the amount of depletant required to induce phase separation 
drops to rjd 8%, showing that adding a small amount of polymer to a binary mixture of hard spheres may give rise 
to sharp demixing provided the polymer-colloid interaction is not too harsh. By further softening the colloid-polymer 
interaction the effect becomes more and more pronounced. A softness exponent n = 12 gives the same sequence of 
phase diagrams but the amount of depletant required to drive the demixing of the colloids is lower than 2%. 

The nature of the transition can be better appreciated by contrasting the shape of the transition lines in the total 
density p = pi + P 2 vs. concentration plane (cj = ^), at a given polymer volume fraction, for two different topologies 
of the phase diagram. As an illustrative example, in Fig. IHlwe compare the case (72 = 0.8 tri with the case (72 = 0.9 (7i 
for the choice ^ = 0.05 cti and n = 36 at 77 ^ = 7.3%. 

As shown in Fig. |S1 although the change in size ratio is rather small (10%), it is enough to switch between the 
two classes of phase diagrams: from type III to type II. In both cases the two-phase region fills a large portion of 
the density-concentration plane, but the nature of the transition changes. To visualize this difference it is instructive 
to follow the tie lines (dotted lines in the figures) which join the two coexisting phases at a given average density p 
and concentration c. It is apparent how the topology of the phase diagram affects the composition of the coexisting 
phases: In type III mixtures (like the case shown in the left panel), when the average total volume fraction is suitably 
chosen, the coexisting phases have both large density and concentration mismatch, while in type I or II systems (right 
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FIG. 8: Same as Fig. [5] for a colloid-depletant interaction with softness exponent n = 36 and ^ = 0.05 cri. The red arrows in 
the left and right panel show the critical points marked by a red dot in Fig. 




FIG. 9: Two slices to the phase diagram in the (p, c) plane for two choices of the colloid size ratio. Here ^ = O.OSci and n = 36. 
The blue line is the spinodal, the full black line is a portion of the coexistence curve and the red dot is the critical point. The 
dotted lines are tie lines, joining the coexisting densities and compositions. The squares (magenta) show the state points where 
correlations have been depicted in Figs. IIOIIII 


panel) the difference in concentration between the coexisting phases is moderate. 

Finally, we briefly discuss how the topology of the phase diagram reflects in particle correlations. The partial 
structure factors are a direct measure of density fluctuations and can be expressed as the equilibrium average of a 
product of density operators of species i and? : S., (k) = ^ (pi(k)pj (— k)) and are related to the direct correlation 
functions by the Ornstein-Zernike relation [J. The structure factors are evaluated by use of the modified RPA 
approximation Eq. da. To discriminate between demixing and liquid-vapour transitions in binary mixtures it is 
convenient to introduce suitable linear combinations of ^^(k), the Bhatia-Thornton structure factors [^ . which 
properly characterize the density and concentration fluctuations: 

Spp = ci5n(k)+C2 522(k)+2(ciC2)^/"5i2(k) 

Sec = cici5ii(k)+C?C2 522(k)-2(ciC2)3/"5l2(k) 

Spe = Sep = Cl C2 [Sll(k) - 522(k)] + (C2 - Cl) (cicz)^/" 5l2(k) (18) 

where Ci = pi/p, p = pi + p 2 are the number concentrations of the two species and the total density, respectively. 

As an illustrative example, the Bhatia-Thornton partial structure factors are shown at the same point in the phase 
diagram for two choices of the size ratio in Figs. (TU] The chosen state point {pa\ = 0.3,ci = 0.1) is marked by 
a square in the sections of the density-concentration plane in Figs. |9l In both cases density correlations appear 
to be significant, even at a rather low volume fraction. Oscillations in Spp{\f.) reveal the presence of short range 
correlations among colloids. At size ratio 0.9, the presence of the liquid vapour critical point (shown in Fig. 12]) 
clearly affects density fluctuations in a large portion of the phase diagram, leading to a strong peak at vanishing wave 
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FIG. 10: Bhatia-Thornton partial structure factors for two choices of the colloid size ratio. The colloid-depletant repulsion is 
chosen to be power-law (HSl): n = 36, ^ = 0.05(71. Density fluctuations (na are shown by the full black curve. Concentration 
fluctuations m by a blue full line, while the dotted red curve gives the off-diagonal term. The density-density structure factor 
for a hard sphere mixture at the same density and concentration is shown by the dashed line. 




FIG. 11: Bhatia-Thornton partial structure factors for two choices of the colloid size ratio. Notation as in Fig. [TO] 


vector. Concentration fluctuations are instead very weak and remarkably flat at this point in the phase diagram. Figs. 
1101 also include the density-density structure factor (1181) for a binary hard sphere mixture at the same density and 
concentration (concentration fluctuations are small and structureless and are not reported). Attractive interactions 
give rise to an enhancement of fluctuations at small wavevectors but also to a shift of the main peak towards higher 
wave vectors, reflecting the tendency to aggregation. These data can be contrasted with the results shown in Fig. 
fm referring to a state point at higher density {paf = 0.6,ci = 0.27), also marked by a square in Figs. HI In both 
cases density correlations acquire a significant structure, showing strong short range correlations induced by depletion, 
particularly at size ratio 0.9 where a very high peak, probably overemphasized by our modified RPA approximation, 
develops. For cr 2 = 0.8 tri, the closeness to the demixing critical point enhances both density and concentration 
fluctuations at long wavelengths. This finding agrees with previous remarks about the character of the demixing 
transition in these systems: the coexisting phases in fact differ both in density and in concentration, implying a 
composite order parameter involving both kinds of fluctuations. In contrast, in a pure hard sphere mixture both 
density and concentration fluctuations are small. 


V. CONCLUSIONS 

The attractive nature of the so-called depletion interaction is well known as well as its effect in the densification of 
colloidal particles. Less well known is that depletion interaction has a demixing effect when colloids of unequal size 



















11 


are present. We have shown how rich is the phase diagram of a mixture of colloids of different but comparable size in 
presence of a depleting agent like small non-adsorbing polymers. In fact the depletion interaction between two spheres 
of different size is sub-additive (see Eq. dH)) compared to the two interactions between pair of spheres of the same 
diameter. Therefore in a binary mixture of colloids two competing effects are present. On one hand the attractive 
nature of the depletion interaction brings in the tendency to phase separation of the colloids in two phases, one rich 
of colloids and one poor of colloids. On the other hand the sub-additivity of the depletion interaction between spheres 
of different diameter promotes phase separation of spheres of different diameter. All three main topologies of the 
phase diagram of a binary mixture described in the Introduction are found, depending on the values of the interaction 
parameters. The tendency to demixing is enhanced by an increasing asymmetry of the colloids and by a decreasing 
size of the polymers. In addition, a strong enhancement of demixing is found if the polymer-colloid interaction has 
some degree of softness and it is not the hard core one assumed in the AO model. In fact, we find that under suitable 
conditions even a concentration of polymers as low as few percent in volume fraction can be enough to cause demixing 
of the colloids. 

We have considered mixtures of colloids of two different sizes. It is clear that our study is relevant also for real 
colloids when the size distribution P{a) is strongly bimodal, i.e. P{cr) as function of the diameter a has two sharp 
maxima. We are not aware of experimental studies of similar colloidal mixtures in the fluid phase in presence of small 
polymers but only at very large volume fraction such colloids has been investigated [2l| with the purpose of studying 
the crystalline phases. In this work a remarkable effect of the polymers on the kinetic of the crystallization was 
found. A binary mixture of hard spheres of similar diameters cti and (72 is fully miscible in the flu id p hase but upon 
solidification phase separation between different crystalline phases can be present . Experiments with mixtures 
of colloids of two sizes at high packing fraction indeed show solidihcation but usually a glass state is formed and the 
kinetic of crystallization is very slow. By adding small polymers in the system one finds a dramatic acceleration 
of the crystallization process and it has been conjectured that this is due to enhanced concentration fluctuations 
induced by the polymer. In that work the size ratio cril<y\ is 0.79 and the polymer diameter ratio ^/(7i = 0.068, 
values similar to those of our study. We can give now a rational basis to the finding of this experiment. Due to the 
sub-additivity of the depletion interaction between spheres of different diameter already in the fluid phase there is 
a tendency to demixing with enhanced concentration fluctuations so that the crystallization in solids with different 
composition should be facilitated. Actually a demixing transition might already be present in the fluid phase. An 
experimental verification of our prediction of demixing in the fluid phase and on the different topology of the phase 
diagram, depending on the interaction parameters, should be very interesting. 

Another interesting issue is how the demixing tendency of the depletion interaction acts at high packing fraction 
of binary colloids where gelation can set in the monodisperse system. 

Our study is based on a rather simple free energy functional of the mixture, it is of a mean field form suitably 
generalized to treat strong and short range attractive forces that are typical of depletion by small particles. Our free 
energy is such that the Noro-Frenkel law of corresponding states is extended to binary mixtures, i.e. thermodynamic 
properties do not depend on the details of the three interaction potentials but only on the three stickiness parameters 
Tij ihj = 1; 2) as given by Eq. ([5]). The theory could be improved but we believe that our basic findings are at least 
qualitatively correct. It should be interesting to check by simulation their accuracy on a quantitative footing. 

Often in the paper we have mentioned non-adsorbing polymers as depleting agent. It is clear that the same 
effects are expected when depletion is induced by other kind of small particles like small micelles. It remains to be 
studied how the sub-additivity of the depletion interaction is modified when the interaction between the depleting 
particles cannot be neglected as it is done in the AO model. An interesting result of our study is the prediction that 
the demixing effect of the depletion interaction is strongly enhanced when some softness is present in the polymer- 
colloid interaction. An experimental study of this aspect should be possible, for instance by coating the colloids by 
adsorbing polymers of different length or by temperature sensitive polymers so that the softness could be controlled 
by temperature. Softening the direct colloid-colloid repulsive interaction may result either in favoring or in hampering 
particle demixing, depending whether the repulsive tail in v‘^'^(r) prevails over the attraction induced by depletion. A 
realistic representation of the pair potential is therefore required to investigate the competition between these effects. 

It is a pleasure to dedicate this work to Jean-Pierre, whose deep insight and outstanding achievements permeate 
modern Liquid State Theory. We acknowledge fruitful discussions with P. Tartaglia, E. Zaccarelli, D. Ashton, N. 
Wilding and A. Giacometti. One of us (L.R.) wants to thank Dipartimento di Fisica, Universita degli Studi di 
Milano, for some support to his research activity. 
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